Comprehensive Transcriptomic Investigation of Rett Syndrome Reveals Increasing Complexity Trends from Induced Pluripotent Stem Cells to Neurons with Implications for Enriched Pathways.

Rett syndrome (RTT) is a rare genetic neurodevelopmental disorder that has no cure apart from symptomatic treatments. While intense research efforts are required to fulfill this unmet need, the fundamental challenge is to obtain sufficient patient data. In this study, we used human transcriptomic data of four different sample types from RTT patients including induced pluripotent stem cells, differentiated neural progenitor cells, differentiated neurons, and postmortem brain tissues with an increasing in vivo-like complexity to unveil specific trends in gene expressions across the samples. Based on DEG analysis, we identified F8A3, CNTN6, RPE65, and COL19A1 to have differential expression levels in three sample types and also observed previously reported genes such as MECP2, FOXG1, CACNA1G, SATB2, GABBR2, MEF2C, KCNJ10, and CUX2 in our study. Considering the significantly enriched pathways for each sample type, we observed a consistent increase in numbers from iPSCs to NEUs where MECP2 displayed profound effects. We also validated our GSEA results by using single-cell RNA-seq data. In WGCNA, we elicited a connection among MECP2, TNRC6A, and HOXA5. Our findings highlight the utility of transcriptomic analyses to determine genes that might lead to therapeutic strategies.

The complete chloroplast genome of Cicer reticulatum and comparative analysis against relative Cicer species.

The chloroplast (cp) genome is an adequate genomic resource to investigate evolutionary relationships among plant species and it carries marker genes available for species identification. The Cicer reticulatum is one of perennial species as the progenitor of cultivated chickpeas. Although a large part of the land plants has a quadruple chloroplast genome organization, the cp genome of C. reticulatum consists of one LSC (Large Single Copy Region), one SSC (Small Single Copy Region), and one IR (Inverted Repeat) region, which indicates that it has an untypical and unique structure. This type of chloroplast genome belongs to the IR-lacking clade. Chloroplast DNA (cpDNA) was extracted from fresh leaves using a high salt-based protocol and sequencing was performed using DNA Nanoball Sequencing technology. The comparative analysis employed between the species to examine genomic differences and gene homology. The study also included codon usage frequency analysis, hotspot divergence analysis, and phylogenetic analysis using various bioinformatics tools. The cp genome of C. reticulatum was found 125,794 bp in length, with an overall GC content of 33.9%. With a total of 79 protein-coding genes, 34 tRNA genes, and 4 rRNA genes. Comparative genomic analysis revealed 99.93% similarity between C. reticulatum and C. arietinum. Phylogenetic analysis further indicated that the closest evolutionary relative to C. arietinum was C. reticulatum, whereas the previously sequenced wild Cicer species displayed slight distinctions across their entire coding regions. Several genomic regions, such as clpP and ycf1, were found to exhibit high nucleotide diversity, suggesting their potential utility as markers for investigating the evolutionary relationships within the Cicer genus. The first complete cp genome sequence of C. reticulatum will provide novel insights for future genetic research on Cicer crops.

Complete chloroplast genome of Lens lamottei reveals intraspecies variation among with Lens culinaris.

Lens lamottei is a member of the Fabaceae family and the second gene pool of the genus Lens. The environmental factors that drove the divergence among wild and cultivated species have been studied extensively. Recent research has focused on genomic signatures associated with various phenotypes with the acceleration of next-generation techniques in molecular profiling. Therefore, in this study, we provide the complete sequence of the chloroplast genome sequence in the wild Lens species L. lamottei with a deep coverage of 713 × next-generation sequencing (NGS) data for the first time. Compared to the cultivated species, Lens culinaris, we identified synonymous, and nonsynonymous changes in the protein-coding regions of the genes ndhB, ndhF, ndhH, petA, rpoA, rpoC2, rps3, and ycf2 in L. lamottei. Phylogenetic analysis of chloroplast genomes of various plants under Leguminosae revealed that L. lamottei and L. culinaris are closest to one another than to other species. The complete chloroplast genome of L. lamottei also allowed us to reanalyze previously published transcriptomic data, which showed high levels of gene expression for ATP-synthase, rubisco, and photosystem genes. Overall, this study provides a deeper insight into the diversity of Lens species and the agricultural importance of these plants through their chloroplast genomes.

Complete chloroplast genome sequence of Lens ervoides and comparison to Lens culinaris.

Lens is a member of the Papilionoideae subfamily of Fabaceae and is generally used as a source of vegetable protein as part of human diets in many regions worldwide. Chloroplast (cp) genomes are highly active genetic components of plants and can be utilized as molecular markers for various purposes. As one of the wild lentil species, the Lens ervoides cp genome has been sequenced for the first time in this study using next-generation sequencing. The de novo assembly of the cp genome resulted in a single 122,722 bp sequence as two separate coexisting structural haplotypes with similar lengths. Results indicated that the cp genome of L. ervoides belongs to the inverted repeat lacking clade. Several noteworthy divergences within the coding regions were observed in ndhB, ndhF, rbcL, rpoC2, and ycf2 genes. Analysis of relative synonymous codon usage showed that certain genes, psbN, psaI, psbI, psbE, psbK, petD, and ndhC, preferred using biased codons more often and therefore might have elevated expression and translation efficiencies. Overall, this study exhibited the divergence level between the wild-type and cultured lentil cp genomes and pointed to certain regions that can be utilized as distinction markers for various goals.

The Complete Chloroplast Genome Sequence of Cicer bijugum, Genome Organization, and Comparison with Related Species.

Background: Chickpea is one of Turkey's most significant legumes, and because of its high nutritional value, it is frequently preferred in human nourishment.Chloroplasts, which have their own genetic material, are organelles responsible for photosynthesis in plant cells and their genome contains non-trivial information about the molecular features and evolutionary process of plants. Objective: Current study aimed at revealing complete chloroplast genome sequence of one of the wild type Cicer species, Cicer bijugum, and comparing its genome with cultivated Cicer species, Cicer arietinum, by using bioinformatics analysis tools. Except for Cicer arietinum, there has been no study on the chloroplast genome sequence of Cicer species.Therefore, we targeted to reveal the complete chloroplast genome sequence of wild type Cicer species, Cicer bijugum, and compare the chloroplast genome of Cicer bijugum with the cultivated one Cicer arietinum. Methods: In this study, we sequenced the whole chloroplast genome of Cicer bijugum, one of the wild types of chickpea species, with the help Next Generation Sequencing platform and compared it with the chloroplast genome of the cultivated chickpea species, Cicer arietinum, by using online bioinformatics analysis tools. Results: We determined the size of the chloroplast genome of C. bijugum as 124,804 bp and found that C. bijugum did not contain an inverted repeat region in its chloroplast genome. Comparative analysis of the C. bijugum chloroplast genome uncovered thirteen hotspot regions (psbA, matK, rpoB, rpoC1, rpoC2, psbI, psbK, accD, rps19, ycf2, ycf1, rps15, and ndhF) and seven of them (matK, accD, rps19, ycf1, ycf2, rps15 and ndhF) could potentially be used as strong molecular markers for species identification. It has been determined that C. bijugum was phylogenetically closer to cultivated chickpea as compared to the other species. Conclusion: It is aimed that the data obtained from this study, which is the first study in which whole chloroplast genomes of wild chickpea species were sequenced, will guide researchers in future molecular, evolutionary, and genetic engineering studies with chickpea species.

HieRFIT: a hierarchical cell type classification tool for projections from complex single-cell atlas datasets.

MOTIVATION: The emergence of single-cell RNA sequencing (scRNA-seq) has led to an explosion in novel methods to study biological variation among individual cells, and to classify cells into functional and biologically meaningful categories. RESULTS: Here, we present a new cell type projection tool, Hierarchical Random Forest for Information Transfer (HieRFIT), based on hierarchical random forests. HieRFIT uses a priori information about cell type relationships to improve classification accuracy, taking as input a hierarchical tree structure representing the class relationships, along with the reference data. We use an ensemble approach combining multiple random forest models, organized in a hierarchical decision tree structure. We show that our hierarchical classification approach improves accuracy and reduces incorrect predictions especially for inter-dataset tasks which reflect real-life applications. We use a scoring scheme that adjusts probability distributions for candidate class labels and resolves uncertainties while avoiding the assignment of cells to incorrect types by labeling cells at internal nodes of the hierarchy when necessary. AVAILABILITY AND IMPLEMENTATION: HieRFIT is implemented as an R package, and it is available at (https://github.com/yasinkaymaz/HieRFIT/releases/tag/v1.0.0). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evaluating single-cell cluster stability using the Jaccard similarity index.

MOTIVATION: One major goal of single-cell RNA sequencing (scRNAseq) experiments is to identify novel cell types. With increasingly large scRNAseq datasets, unsupervised clustering methods can now produce detailed catalogues of transcriptionally distinct groups of cells in a sample. However, the interpretation of these clusters is challenging for both technical and biological reasons. Popular clustering algorithms are sensitive to parameter choices, and can produce different clustering solutions with even small changes in the number of principal components used, the k nearest neighbor and the resolution parameters, among others. RESULTS: Here, we present a set of tools to evaluate cluster stability by subsampling, which can guide parameter choice and aid in biological interpretation. The R package scclusteval and the accompanying Snakemake workflow implement all steps of the pipeline: subsampling the cells, repeating the clustering with Seurat and estimation of cluster stability using the Jaccard similarity index and providing rich visualizations. AVAILABILITYAND IMPLEMENTATION: R package scclusteval: https://github.com/crazyhottommy/scclusteval Snakemake workflow: https://github.com/crazyhottommy/pyflow_seuratv3_parameter Tutorial: https://crazyhottommy.github.io/EvaluateSingleCellClustering/.

Epstein-Barr Virus Genomes Reveal Population Structure and Type 1 Association with Endemic Burkitt Lymphoma.

Endemic Burkitt lymphoma (eBL), the most prevalent pediatric cancer in sub-Saharan Africa, is distinguished by its inclusion of Epstein-Barr virus (EBV). In order to better understand the impact of EBV variation in eBL tumorigenesis, we improved viral DNA enrichment methods and generated a total of 98 new EBV genomes from both eBL cases (n = 58) and healthy controls (n = 40) residing in the same geographic region in Kenya. Using our unbiased methods, we found that EBV type 1 was significantly more prevalent in eBL patients (74.5%) than in healthy children (47.5%) (odds ratio = 3.24, 95% confidence interval = 1.36 to 7.71, P = 0.007), as opposed to similar proportions in both groups. Controlling for EBV type, we also performed a genome-wide association study identifying six nonsynonymous variants in the genes EBNA1, EBNA2, BcLF1, and BARF1 that were enriched in eBL patients. In addition, viruses isolated from plasma of eBL patients were identical to their tumor counterparts consistent with circulating viral DNA originating from the tumor. We also detected three intertypic recombinants carrying type 1 EBNA2 and type 2 EBNA3 regions, as well as one novel genome with a 20-kb deletion, resulting in the loss of multiple lytic and virion genes. Comparing EBV types, viral genes displayed differential variation rates as type 1 appeared to be more divergent, while type 2 demonstrated novel substructures. Overall, our findings highlight the complexities of the EBV population structure and provide new insight into viral variation, potentially deepening our understanding of eBL oncogenesis.IMPORTANCE Improved viral enrichment methods conclusively demonstrate EBV type 1 to be more prevalent in eBL patients than in geographically matched healthy controls, which previously underrepresented the prevalence of EBV type 2. Genome-wide association analysis between cases and controls identifies six eBL-associated nonsynonymous variants in EBNA1, EBNA2, BcLF1, and BARF1 genes. Analysis of population structure reveals that EBV type 2 exists as two genomic subgroups and was more commonly found in female than in male eBL patients.

Evolutionarily conserved pachytene piRNA loci are highly divergent among modern humans.

In the fetal mouse testis, PIWI-interacting RNAs (piRNAs) guide PIWI proteins to silence transposons but, after birth, most post-pubertal pachytene piRNAs map to the genome uniquely and are thought to regulate genes required for male fertility. In the human male, the developmental classes, precise genomic origins and transcriptional regulation of postnatal piRNAs remain undefined. Here, we demarcate the genes and transcripts that produce postnatal piRNAs in human juvenile and adult testes. As in the mouse, human A-MYB drives transcription of both pachytene piRNA precursor transcripts and messenger RNAs encoding piRNA biogenesis factors. Although human piRNA genes are syntenic to those in other placental mammals, their sequences are poorly conserved. In fact, pachytene piRNA loci are rapidly diverging even among modern humans. Our findings suggest that, during mammalian evolution, pachytene piRNA genes are under few selective constraints. We speculate that pachytene piRNA diversity may provide a hitherto unrecognized driver of reproductive isolation.

Integrative microRNA and mRNA deep-sequencing expression profiling in endemic Burkitt lymphoma.

BACKGROUND: Burkitt lymphoma (BL) is characterized by overexpression of the c-myc oncogene, which in the vast majority of cases is a consequence of an IGH/MYC translocation. While myc is the seminal event, BL is a complex amalgam of genetic and epigenetic changes causing dysregulation of both coding and non-coding transcripts. Emerging evidence suggest that abnormal modulation of mRNA transcription via miRNAs might be a significant factor in lymphomagenesis. However, the alterations in these miRNAs and their correlations to their putative mRNA targets have not been extensively studied relative to normal germinal center (GC) B cells. METHODS: Using more sensitive and specific transcriptome deep sequencing, we compared previously published small miRNA and long mRNA of a set of GC B cells and eBL tumors. MiRWalk2.0 was used to identify the validated target genes for the deregulated miRNAs, which would be important for understanding the regulatory networks associated with eBL development. RESULTS: We found 211 differentially expressed (DE) genes (79 upregulated and 132 downregulated) and 49 DE miRNAs (22 up-regulated and 27 down-regulated). Gene Set enrichment analysis identified the enrichment of a set of MYC regulated genes. Network propagation-based method and correlated miRNA-mRNA expression analysis identified dysregulated miRNAs, including miR-17~95 cluster members and their target genes, which have diverse oncogenic properties to be critical to eBL lymphomagenesis. Central to all these findings, we observed the downregulation of ATM and NLK genes, which represent important regulators in response to DNA damage in eBL tumor cells. These tumor suppressors were targeted by multiple upregulated miRNAs (miR-19b-3p, miR-26a-5p, miR-30b-5p, miR-92a-5p and miR-27b-3p) which could account for their aberrant expression in eBL. CONCLUSION: Combined loss of p53 induction and function due to miRNA-mediated regulation of ATM and NLK, together with the upregulation of TFAP4, may be a central role for human miRNAs in eBL oncogenesis. This facilitates survival of eBL tumor cells with the IGH/MYC chromosomal translocation and promotes MYC-induced cell cycle progression, initiating eBL lymphomagenesis. This characterization of miRNA-mRNA interactions in eBL relative to GC B cells provides new insights on miRNA-mediated transcript regulation in eBL, which are potentially useful for new improved therapeutic strategies.

FLT1 and transcriptome-wide polyadenylation site (PAS) analysis in preeclampsia.

Maternal symptoms of preeclampsia (PE) are primarily driven by excess anti-angiogenic factors originating from the placenta. Chief among these are soluble Flt1 proteins (sFlt1s) produced from alternatively polyadenylated mRNA isoforms. Here we used polyadenylation site sequencing (PAS-Seq) of RNA from normal and PE human placentae to interrogate transcriptome-wide gene expression and alternative polyadenylation signatures associated with early-onset PE (EO-PE; symptom onset < 34 weeks) and late-onset PE (LO-PE; symptom onset > 34 weeks) cohorts. While we observed no general shift in alternative polyadenylation associated with PE, the EO-PE and LO-PE cohorts do exhibit gene expression profiles distinct from both each other and from normal placentae. The only two genes upregulated across all transcriptome-wide PE analyses to date (microarray, RNA-Seq and PAS-Seq) are NRIP1 (RIP140), a transcriptional co-regulator linked to metabolic syndromes associated with obesity, and Flt1. Consistent with sFlt1 overproduction being a significant driver of clinical symptoms, placental Flt1 mRNA levels strongly correlate with maternal blood pressure. For Flt1, just three mRNA isoforms account for > 94% of all transcripts, with increased transcription of the entire locus driving Flt1 upregulation in both EO-PE and LO-PE. These three isoforms thus represent potential targets for therapeutic RNA interference (RNAi) in both early and late presentations.

High pathogen burden in childhood promotes the development of unconventional innate-like CD8+ T cells.

Cellular and humoral constituents of the immune system differ significantly between children and adults, yet very little is known about the impact of early-life pathogen exposure on this immunologic transition. We examined CD4+ and CD8+ T cell subsets defined by CCR7 and CD45RA expression in two longitudinal pediatric cohorts experiencing divergent levels of pathogen burden. Using multiparameter flow cytometry, along with serological, cytokine, and transcriptomic data, we show that cumulative pathogen burden promotes the development of atypical CD8dim T cells with an innate-like profile (Granzyme Bhi, IFNγlow, TNFαlow, PLFZhi, ID2hi, IKZF2hi) in contrast to age-matched children residing in a low pathogen-exposure area who display a more conventional CD8bright profile (IFNγ+, TNFα+, CCL4+). Furthermore, these unconventional T cells had stunted proliferation, distinct transcriptional programs, and impaired T cell receptor signaling and were enriched in hallmark TNFα, NF-κB, and IL-6 gene signaling pathways, reminiscent of NK cells and type-1 innate lymphoid cells. Our findings suggest that these unconventional CD8dim T cells arise in a very particular immunological context and may provide a deeper understanding of the heterogeneity in human immune responses.

Comprehensive Transcriptome and Mutational Profiling of Endemic Burkitt Lymphoma Reveals EBV Type-Specific Differences.

Endemic Burkitt lymphoma (eBL) is the most common pediatric cancer in malaria-endemic equatorial Africa and nearly always contains Epstein-Barr virus (EBV), unlike sporadic Burkitt lymphoma (sBL) that occurs with a lower incidence in developed countries. Given these differences and the variable clinical presentation and outcomes, we sought to further understand pathogenesis by investigating transcriptomes using RNA sequencing (RNAseq) from multiple primary eBL tumors compared with sBL tumors. Within eBL tumors, minimal expression differences were found based on: anatomical presentation site, in-hospital survival rates, and EBV genome type, suggesting that eBL tumors are homogeneous without marked subtypes. The outstanding difference detected using surrogate variable analysis was the significantly decreased expression of key genes in the immunoproteasome complex (PSMB9/ß1i, PSMB10/ß2i, PSMB8/ß5i, and PSME2/PA28ß) in eBL tumors carrying type 2 EBV compared with type 1 EBV. Second, in comparison with previously published pediatric sBL specimens, the majority of the expression and pathway differences was related to the PTEN/PI3K/mTOR signaling pathway and was correlated most strongly with EBV status rather than geographic designation. Third, common mutations were observed significantly less frequently in eBL tumors harboring EBV type 1, with mutation frequencies similar between tumors with EBV type 2 and without EBV. In addition to the previously reported genes, a set of new genes mutated in BL, including TFAP4, MSH6, PRRC2C, BCL7A, FOXO1, PLCG2, PRKDC, RAD50, and RPRD2, were identified. Overall, these data establish that EBV, particularly EBV type 1, supports BL oncogenesis, alleviating the need for certain driver mutations in the human genome. IMPLICATIONS: Genomic and mutational analyses of Burkitt lymphoma tumors identify key differences based on viral content and clinical outcomes suggesting new avenues for the development of prognostic molecular biomarkers and therapeutic interventions.

Human and Epstein-Barr Virus miRNA Profiling as Predictive Biomarkers for Endemic Burkitt Lymphoma.

Endemic Burkitt lymphoma (eBL) is an aggressive B cell lymphoma and is associated with Epstein-Barr virus (EBV) and Plasmodium falciparum malaria co-infections. Central to BL oncogenesis is the over-expression of the MYC proto-oncogene which is caused by a translocation of an Ig enhancer in approximation to the myc gene. While whole genome/transcriptome sequencing methods have been used to define driver mutations and transcriptional dysregulation, microRNA (miRNA) dysregulation and differential expression has yet to be fully characterized. We hypothesized that both human and EBV miRNAs contribute to eBL clinical presentation, disease progression, and poor outcomes. Using sensitive and precise deep sequencing, we identified miRNAs from 17 Kenyan eBL patient tumor samples and delineated the complement of both host and EBV miRNAs. One human miRNA, hsa-miR-10a-5p was found to be differentially expressed (DE), being down-regulated in jaw tumors relative to abdominal and in non-survivors compared to survivors. We also examined EBV miRNAs, which made up 2.7% of the miRNA composition in the eBL samples. However, we did not find any significant associations regarding initial patient outcome or anatomical presentation. Gene ontology analysis and pathway enrichment of previously validated targets of miR-10a-5p suggest that it can promote tumor cell survival as well as aid in evasion of apoptosis. To examine miR-10a-5p regulatory effect on gene expression in eBL, we performed a pairwise correlation coefficient analysis on the expression levels of all its validated targets. We found a significant enrichment of correlated target genes consistent with miR-10a-5p impacting expression. The functions of genes and their correlation fit with multiple target genes impacting tumor resilience. The observed downregulation of miR-10a and associated genes suggests a role for miRNA in eBL patient outcomes and has potential as a predictive biomarker that warrants further investigation.

Biodegradation kinetics of o-cresol by Pseudomonas putida DSM 548 (pJP4) and o-cresol removal in a batch-recirculation bioreactor system

The biodegradation kinetics of o-cresol was examined by acclimatized P. putida DSM 548 (pJP4) in batch experiments at varying initial o-cresol concentrations (from 50 to 500 mg/L). The kinetic parameters of o-cresol aerobic biodegradation were estimated by using the Haldane substrate inhibition equation. The biodegradation kinetics of o-cresol was investigated. In batch culture reactors, the Maximum specific growth rate (μmax), Monod constant (Ks) and the inhibition constant (Ki) were established as 0.519 h-1, 223.84 mg/L and 130.883 mg/L, respectively. o-cresol biodegradation in a batch-recirculation bioreactor system by immobilized P. putida was also studied. The recycled packed bed reactor system, which was composed of Ca-alginate beads and pumice on which cells immobilized, has been performed to determine possible stability for further developments.

Biosorption kinetics of a direct azo dye Sirius Blue K-CFN by Trametes versicolor

In this study, lyophilized Trametes versicolor biomass is used as a sorbent for biosorption of a textile dye, Sirius Blue K-CFN, from an aqueous solution. The batch sorption was studied with respect to dye concentration, adsorbent dose and equilibrium time. The effect of pH and temperature on dye uptake was also investigated and kinetic parameters were determined. Optimal initial pH (3.0), equilibrium time (2 hrs), initial dye concentration ( 100 mg l-1) and biomass concentration (1.2 mg l-1) were determined at 26ºC. The maximum biosorption capacity (q max) of Sirius Blue K-CFN dye on lyophilized T. versicolor biomass is 62.62 mg/g. The kinetic and isotherm studies indicated that the biosorption process obeys to a pseudo-second order model and Langmuir isotherm model. In addition, the biosorption capacities of fungal biomass compared to other well known adsorbents such as activated carbon and Amberlite, fungal biomass biosorptions capacities were found to be more efficient.